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A high-resolution calorimetric study has been carried out on nano-coUoidal dispersions of aerosils 
in the liquid crystal 4-n-pentylphenylthiol-4'-n-octyloxybenzoate (8S5) as a function of aerosil con- 
centration and temperature spanning the smectic-C to nematic phases. Over this temperature 
range, this liquid crystal possesses two continuous XY phase transitions: a fluctuation dominated 
nematic to smectic-^ transition with a ~ Oxy ~ —0.013 and a mean-field smectic-^ to smectic-C 
transition. The effective critical character of the A'^-Smj4 transition remains unchanged over the 
entire range of introduced quenched random disorder while the peak height and enthalpy can be 
well described by considering a cut-off length scale to the quasi-critical fluctuations. The robust 
nature of the A''-Sm^ transition in this system contrasts with cyanobiphenyl-aerosil systems and 
may be due to the mesogens being non-polar and having a long nematic range. The character of 
the SmA-SmC transition changes gradually with increasing disorder but remains mean-field- like. 
The heat capacity maximum at the Smj4-SmC transition scales as p^"'^ with an apparent evolution 
from tricritical to a simple mean-field step behavior. These results may be generally understood as 
a stiffening of the liquid crystal (both the nematic elasticity as well as the smectic layer compression 
modulus B) with silica density. 

PACS numbers: 64.70.Md, 61.30.Eb, 65.40.Ba 



INTRODUCTION 



The study of quenched random disorder effects ad- 
dresses many fundamental issues of current interest in 
statistical mechanics. Recent experimental advances 
have shed considerable light onto the random-field theo- 
retic approach, believed to be underlying the physics of 
quenched random disorder (QRD). In particular, these 
efforts have led to the systematic study of the random- 
field model for transitions that break a continuous sym- 
metry 01 . This model is applicable to a terrific range of 
phenomena: This include unique assemblies of colloids, 
complex fluids, charge density waves, spin glasses, and 
doped semiconductors, to name a few. The experimental 
efforts to date have focussed on idealized physical systems 
in order to isolate the essential features of quenched ran- 
dom disorder. Considerable research has been carried out 
on the superfluid transition of ^He and ^He-^He mixtures 
in a variety of porous media 0, 13 as well as doped mag- 
netic systems |J| but many questions remain due to the 
quantum nature of the former and the glass-like behavior 
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of the latter 0. 

In the random- field approach, the effect of the disor- 
der is mapped onto a local field, h{f), directly coupled to 
the order parameter. This local field varies randomly 
through the system on length scales smaller than the 
length scales of the ordered phase such that {h)f — 
while {h ■ h) f ^ 0. Because the nature of the imposed 
disorder is modelled as a dilute (or weak) random-field, 
the effect of QRD may be understood statistically. By 
studying in detail good realizations of particular univer- 
sality classes, the results can be applied more generally to 
a wide variety of physical systems. A particularly fruitful 
physical system to explore QRD effects in general and 
the random-field model in particular, have been nano- 
coUoidal dispersions of an aerosil gel in a liquid crystal 
(LC) host (LC-hsil) [l^. 

The most well studied LC+sil phase transition has 
been the continuous nematic to smectic-j4 (A^-Smj4) 
phase transition. This transition involves the break- 
ing of a continuous symmetry and belongs, though not 
simply, to the 3D-XY universality class In gen- 

eral, high-resolution x-ray 0, H, 0, 0| and calorime- 
try ^.lli 12, 13, 14, ISj studies on the iV-Sm^+sil tran- 
sition have found that the (quasi-) long-range ordered 
smectic phase is destroyed for all densities of aerosil. 



2 



However, below a silica density of about 0.1 grams of 
SiOq per cni^ of LC (the so-called conjugate density, 
PS whose units will be dropped hereafter), a pseudo- 
transition displaying quasi-critical behavior ai T — T* 
persists. The scattering from the short-range smectic 
order is well described by both a bulk thermal struc- 
ture factor (dominate above the pseudo-transition T*) 
and a random- field structure factor (dominate below T*) 
given by the bulk thermal form squared, the so-called 
" lorentzian-|-lorentzian squared" form 0, 0. Although 
the smectic correlation length is finite, it is large (span- 
ning many mean-void lengths of the silica gel) and ex- 
hibits a power-law divergence on cooling to T* at which 
it saturates for 8CB and 80CB in aerosils or begins 
to decrease upon further cooling towards the SmC phase 
for 8S5-hsil system This broadening observed for 

8S5-hsil [13 and anisotropic (field aligned) 8S5-hsil [l^ 
samples, could be due to a distribution of SmC tilt 
angles and result in a lower average order parameter 
squared. The calorimetry studies found sharp, quasi- 
critical, power-law divergences for the heat capacity for 
ps < 0.1, while all phase transitions appear to be elas- 
tically smeared for larger ps |3. The evolution of an ef- 
fective critical heat capacity exponent cteff with aerosil 
density is consistent with a gradual drop in the nematic 
susceptibility with increasing gel density for all A'^-Smj4 
transitions studied to date. 

A detailed scaling analysis of the A^-Sm^ transition 
combining x-ray and calorimetry results through the 
quasi-critical transition for low aerosil density revealed 
the importance of both random-field and finite-size-like 
effects 1^ 01 . The dimensionality of the system also 
played a role and was consistent with the expected di- 
mensional rise, cL/ = d + 2, for a system with random- 
field disorder |3 ^M- This work also found in gen- 
eral that the effective random-field strength (or variance) 
scales as ps over the entire range of silica densities stud- 
ied thus far, although there are indications of different 
scalings above and below ps = 0.1 |l7j . The experimen- 
tal and theoretical results to date support the view that 
a random-field- XY (RF-XY) system has no new critical 
point. 

Recently, tilted smectic phases have become the focus 
of studies on quenched random disorder effects 0, 0| . 
The smectic-^ to smectic- C (Sm^-SmC) phase transi- 
tion involves the breaking of a continuous symmetry and 
is described by two parameters, the tilt and azimuthal an- 
gles. This transition belongs to the 3D-XY universality 
class but is mean-field in character due to the relatively 
long bare correlation length of the SmC. The strong cou- 
pling between tilt and smectic layer compression for the 
SmC phase appears to place the Sm^-SmC transition al- 
ways close to a (Landau) mean- field tricritical point 6J. 
A consequence of the tilt angle's sensitivity to the layer 
elasticity is that this transition is much more strongly 
disordered by perturbations that distort the smectic lay- 
ers, as in the strong disorder of aerogels, than the A^- 
Sm^ pol l2l|. For weaker aerosil gel induced disorder. 



the Landau mean-field heat capacity signature was found 
to be unaffected by hydrophobic aerosils [Ij] while a re- 
cent high-resolution x-ray scattering study on a different 
LC in a hydrophilic aerosil found little change in the tilt 
angle temperature dependence with ps 10\. This robust- 
ness to weak disorder may be a consequence of the mean- 
field character of the transition placing it effectively at its 
upper critical dimension. However, when the transition 
from the Sm^ phase goes to the chiral analog of the SmC 
phase (SmC*), the effect of even the aerosil induced dis- 
order generally suppresses and smears the transition j^- 

The present work focusses on the effect of quenched 
random disorder induced by a nano-coUoidal dispersion 
of hydrophilic typc-300 aerosil forming a mass-fractal 
gel within the liquid crystal 4-n-pentylphcnylthiol-4'- 
Ti-octyloxybenzoate (8S5) that closely follows the pre- 
viously reported x-ray studies ^Q. High-resolution 
calorimetry has been carried out on these 8S5-|-sil dis- 
persions as a function of aerosil concentration and tem- 
perature spanning the smectic- C to nematic phases. This 
liquid crystal possesses two continuous XY phase tran- 
sitions of interest. The first is a fluctuation dominated 
nematic to smectic-^ characterized by a heat capacity 
exponent a < 0. The critical character of the 7V-Sm^ 
transition remains unchanged with the introduction of 
the quenched random disorder while its enthalpy and 
heat capacity maximum decreases in a manner consis- 
tent with a finite-size-like scaling without any obvious 
crossover from soft to stiff gel seen on other LC-|-sil sys- 
tems. The stability of the heat capacity quasi-critical 
behavior for this system with QRD is consistent with 
our understanding of the influence of the gel on the un- 
derlying nematic order. Since in the bulk material the 
nematic susceptibility is low before Sm^ order begins to 
form, the gel has little scope to reduce this quantity. 

The second is a mean-field Landau tricritical smectic-^ 
to smectic- C phase transition. The Smj4-SmC remains 
mean-field for all aerosil concentrations studied with an 
continuous evolution from the bulk's tricritical to a sim- 
ple mean-field step heat capacity behavior analogous to 
that seen for the pure Sm^-SmC* transition in binary 
LC mixtures In particular, the Sm^l-SmC heat ca- 
pacity maximum at the transition scales as Pg"'^. The 
stable mean-field character of the Sm^-SmC with QRD 
may be a consequence of this transition being effectively, 
due to the long range interactions, at its upper critical 
dimension. The observed crossover from tricritical to a 
simple mean-field step behavior for the Sm^-SmC and 
the continuous diminishing of the A^-Sm^ heat capacity 
peaks may be understood as a continuous stiffening of 
the liquid crystal elasticity for the nematic and smectic 
structure with increasing silica density. However, many 
theoretical challenges are evident. 

Section ^1 describes the preparation of the 8S5+sil 
dispersions as well as the ac-calorimetry technique em- 
ployed. Section IIhI presents the results of the calorimet- 
ric study, while Section IfVI discusses the significance of 
the evolution of an XY transition as well as a crossover 
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from Landau tricritical to mean-field continuous transi- 
tion with increasing quenched disorder. Directions for 
future study will also be discussed. 



II. SAMPLES AND CALORIMETRY 

The liquid crystal 8S5, synthesized at Kent State Uni- 
versity, was used after degassing in the isotropic phase for 
two hours. The best literature reported transition tem- 
perature values in the bulk for this liquid-crystal molecule 
( = 412.64 g mol-i ) are ^ 359.6 K for the 
weakly first-order isotropic to nematic (I-N) transition, 
= 336.58 K for the XY-like continuous nematic to 
smectic-^ (#-Sm^) transition, and T^c - 329.35 K for 
the monotropic Landau tricritical smectic-^ to smectic- 
C (Smj4-SmC) transition 23]. At lower temperatures 
on cooling, a monotropic smectic-C to crystal-i3 (SmC- 
GrB) transition occurs at Tq^ ^ 304 K. The strongly 
first-order Crystal-Smyl (Cr-Sm^) transition occurs re- 
producibly on heating at T^^A ~ 332 K. The mea- 
sured transitions temperatures for our bulk material oc- 
cur at T^^ = 336.64 K, T^^ ^ 328.96 K, and Tg^^ ~ 
331.0 K, which are in reasonable agreement with the 
literature bulk values. However, 8S5 is known to age 
with its transition temperatures continuously shift down- 
ward with time, especially when heated into the isotropic 
phase 23, 24J. Fortunately, the iV-Smvl and Sm^-SmC 
transitions remain relatively sharp, well defined, and con- 
sistent in shape during the sample aging, see ^g.^a-nd 
12 as well as Table U Although the aging of this LC is 
unavoidable, the samples studied in this work all expe- 
rienced the same preparation method and thermal his- 
tory, thus the relative evolution of these transitions with 
aerosil disorder should be preserved. 

The nano-colloidal mixture of hydrophilic type-300 
aerosil in 8S5 was prepared following the solvent disper- 
sion procedure outlined in Refs. P,u^. The hydrophilic 
nature of the aerosils allows the silica particles to weakly 
hydrogen bond to each other and form a gel in an organic 
solvent. The specific surface area of type-300 aerosil is 
300 m^ g~^ [23 and each aerosil sphere is roughly 7 nm 
in diameter. However, the basic free-floating aerosil unit 
consists of a few of these spheres fused together during 
the manufacturing process [lj| • Each 8S5-|-sil sample was 
created by mixing appropriate quantities of liquid-crystal 
and aerosil together, then dissolving the resulting mix- 
ture in spectroscopic grade (low water content) acetone. 
The resulting solution was then dispersed using an ul- 
trasonic bath for about an hour. As the acetone evapo- 
rates from the mixture, a fractal-like gel forms through 
diffusion-limited aggregation. Crystallization of the LC 
host can severely disrupt the gel structure and so care 
was taken to prevent any formation of the solid phase of 
the liquid crystal during the experiments. 

High-resolution ac calorimetry was performed using 
two home-built calorimeters at WPI. The sample cell 
consisted of a silver crimped-sealed envelope ~ 10 mm 




FIG. 1: The specific heat of bulk 8S5 on cooling spanning 
the nematic, smectic-A, and smectic-C phases. The dashed 
line represents the linear Cp (background) used to extract the 
excess Cp associated with the N-SmA transition, ACp{N A) . 
The solid line represents the low-temperature wing of the N- 
SmA transition and is used to extract the excess Cp associated 
with the Smj4-SmC transition, SCp{AC). 
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FIG. 2: Excess specific heat ACp of the A''-Smyl transition 
±3.5 K about Tna for bulk 8S5 and 8S5+siI samples. Inset 
lists the ps for each data set shown. The temperature range 
shown corresponds approximately to ±10"^ width in reduced 
temperature. 
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TABLE I: Summary of the calorimetric results for 8S5+sil samples. Shown are the conjugate silica density (ps in grams of 
aerosil per cm^ of 8S5), the mean- void length lo = 2/aps (where a = 300 m^ g"'^ is the specific surface area of this aerosil) 
within the gel in A, the N-SmA (Tjva) and the Smj4-SmC (Tag) pseudo-phase transition temperatures, and the smectic-A 
temperature range (ATa = Tna — Tac) all in kelvins and averaged between heating and cooling scans. These are followed by 
similarly averaged enthalpy (SHna) in J g~^ and specific heat maximum hm = AC^"^ in J g~^ values for the 7V-Sm^ 
pseudo-phase transition. The final column tabulates the specific heat step in J g~^ of the SmA-SmC phase transition 
^C'AC^ averaged between heating and cooling scans and taken as the value of the excess specific heat —6 K below Tac- 
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long, ^ 5 mm wide, and ~ 0.5 mm thick (closely match- 
ing the dimensions of the heater). After the sample was 
introduced into a cell having an attached 120-il strain- 
gauge heater and l-MfJ carbon-flake thermistor, a con- 
stant current was placed across the heater to maintain 
the cell temperature well above Tjjy. The filled cell was 
then placed in an ultrasonic bath to remix the sample. 
After remixing, the cell was mounted in the calorime- 
ter, the details of which have been described elsewhere 
|2g . In the ac-mode, power is input to the cell as Poce*"* 
resulting in temperature oscillations with amplitude Tac 
and a relative phase shift oi (p = $ + 7r/2, where $ is 
the absolute phase shift between Tad^^) and the input 
power. The specific heat at a heating frequency w using 
C* = Pac/uj \Tac\ is given by 

„ _ filled ~ ^empty] _ C* COSif — C empty /-.x 

P — — — ' Uj 

fl^sample l^sample 



Cfuied = C* sirup - —— , (2) 

where Cj^^^^ and Cj^^^j^ are the real and imaginary com- 
ponents of the heat capacity, Cempty is the heat capacity 
of the cell and silica, nisampie is the mass in grams of the 
liquid crystal (the total mass of the 8S5-|-sil sample was 
~ 20 mg, which yielded msampie values in the range of 
13 — 20 mg), and Re is the external thermal resistance 
linking the cell and the bath (here, ~ 200 K W~-^). Equa- 
tions and lj2Jl require a small correction to account for 
the finite internal thermal resistance compared to R, and 
this was applied to all samples studied here [23|. Mea- 
surements were conducted at various frequencies from 1 
to 500 mHz in order to ensure the applicability of Eqs. ^ 
and (dJ. For all results presented here, Coined — 
through the A^-Sm^ and Smj4-SmC transitions, which 
is expected for continuous phase transitions, and that Cp 
was independent of lu. All data presented here were taken 



at a heating frequency oi ut = 0.1473 s^^ (or 23.4 mHz) 
and a scanning rate of less than ±100 mK h^^, which 
yielded static Cp results. This equilibrium behavior re- 
lates only to the pseudo-critical smectic fluctuations that 
are the focus of this work. The glass-like behavior, that 
has been observed in rheological studies in this frequency 
range |35j | , is a characteristic of these systems away from 
the pseudo-transition." 

The bulk 8S5 and all 8S5-f sil samples experienced the 
same thermal history after mounting; six hours in the 
isotropic phase to ensure homogeneous gelation, then a 
slow cool into the smectic phase at ^ 320 K before begin- 
ning the first detailed heating scan to ^ 345 K followed 
immediately by a detailed cooling scan over the same 
range. Tests on bulk 8S5 with various thermal histories 
in the isotropic phase reveal an aging of the material. 
There are progressive shifts of the transition tempera- 
tures downward with increased time at high tempera- 
ture. However, other than an increase in the rounding 
of the Cp peaks and the downward shift of the transi- 
tion with time, the critical character remains essentially 
unchanged. 

III. RESULTS 

The specific heat on cooling for a bulk sample of 8S5 is 
shown in Fig. ^ Clearly visible are the XY-like A-Smj4 
phase transition at 336.71 K and a Landau mean-field 
Smj4-SmC at 328.82 K. These are in good agreement 
with Schantz and Johnson The excess specific heat 
due to the A-Sm^ transition ACp{NA) for the bulk and 
LC-|-sil samples is obtained by subtracting from the spe- 
cific heat Cp a linear background 

ACp = Cp — Cp (background) , (3) 

where Cp (background) is shown as the dashed line in 
Fig. n and represents the Cp variation of the low- 
temperature wing of the I-N transition. This expression 
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Ar=r-r^^(K) 

FIG. 3; Excess specific heat SCp of tfie Smyl-SmC phase 
transition from +2.5 K above to —4.5 K below Tac for bulk 
8S5 and 8S5+sil samples. Inset lists the ps for each data set 
shown. 

is a valid representation of ACp{NA) for all T > Tac 
and the resulting ACp{NA) for all samples studied are 
shown in Fig. [21 over a ±3.5 K temperature range about 
Tna- The fluctuation dominated A'^-Sm^ transition en- 
thalpy is defined as 

SHna^ J ACpiNA)dT , (4) 

where consistent limits of the integration of ±5 K about 
Tna were used for all samples. The specific heat con- 
tribution of the Sm^-SmC transition is obtained from 
ACp by subtracting the low-temperature heat capacity 
wing of the A''-Sm^ transition (see Fig.nj: 

SCp{AC) = ACp - AC7"9(iVA) . (5) 

The result for all samples studied are shown in Fig. O 
over a range +2 K above to —4.5 K below Tac- The rel- 
evant thermal characteristics for both phase transitions 
of 8S5-|-sils as a function of ps are given in Table 

As seen in Fig. |21 the A^-SmA Cp peak remains sharp 
for all aerosil densities up to the maximum studied of 
PS = 0.347 with no abrupt truncation marking a transi- 
tion from soft to stiff gel behavior seen in 8CB-|-sil 0| 
and 80CB-hsil systems. The ACp{NA) wings on 
both sides of the transition decrease with increasing ps 
in contrast with that found for the 8CB and 80CB in 
aerosils, which found significant reduction of the high- 
temperature ACp{NA) wing Qualitatively, the 
ACp{NA) peak shape remains remarkably stable with 
increasing ps- For the 8S5-|-sil system, the most striking 



effect of the aerosils is the strong decrease in the heat ca- 
pacity maximum at the transition, hm = AC™"^ , 
and in the transition enthalpy SH]\fA with increasing ps- 
In addition, the value of Tna{ps) a-lso displays very lit- 
tle sensitivity to the aerosils even when the aging of the 
8S5 is taken into account, in stark contrast to the non- 
monotonic downward shifts found for nearly all other 
LC+sil systems JJJ. See Table |l| 

The effect of increasing aerosil concentration on the 
Sm^-SmC transition heat capacity appears entirely at 
and below the transition Tac- The 5Cp{AC) peak 
smoothly evolves with increasing ps from the bulk tri- 
critical behavior to a simple heat capacity step as shown 
in Fig.|31 Also apparent is the systematic reduction of the 
step magnitude with increasing P5 as listed in Table ^ 
While there is some rounding of 6Cp{AC) in the imme- 
diate vicinity of Tac, there does not appear any strong 
effect of the aerosil on the high-temperature SCp tail. As 
with the A'^-Sm^ transition, even when aging of the bulk 
material is taken into account, Tac{ps) is not shifted 
downward with increasing aerosil content but appears to 
shift slightly upward. The stability of both Tna{ps) and 
Tac{ps) reported here is in contrast to that found in 
the recent x-ray study that found large non-monotonic 
transition temperature shifts with ps However, the 
thermal history experienced by these calorimetry sam- 
ples was far less severe than that experienced by the x- 
ray samples which were held at high temperatures while 
the solvent evaporated prior to the measurements. Thus, 
the results presented here should reflect a consistent mea- 
sure of the transition temperature shift with ps for this 
system. 

A. The XY-like JV-Sm.4. scaling analysis 

The shape of these experimental ACp{NA) data as a 
function of aerosil content is characterized by a power-law 
form pol l in terms of the reduced temperature t =\ T — 
T* I / T* given by 

ACp{NA)^A^t-°'{\ + D^t^^) + B^, (6) 

where the critical behavior as a function of reduced tem- 
perature t is characterized by an exponent a, amplitudes 
A^ above and below the transition, a critical background 
term and corrections-to-scaling terms having an am- 
plitude and exponent Ai ~ 0.5. An increasing tem- 
perature gap of excluded data about the ACp{NA) peak 
with increasing ps was required to perform the nonlinear 
least-squares fitting. These excluded data were identified 
by strong deviations in the residuals near T* and repre- 
sent rounded ACp{NA) values not describable by the 
power-law given in Eq. © . These fit results for the sets 
of heating and cooling scans of bulk 8S5 and 8S5-|-sil sam- 
ples are presented in Tables ^ and IIIII respectively. In 
addition to fits with the effective critical exponent ae// 
as a free parameter, fits were also performed with Ofe// 
fixed to the 3D-XY value of axY = -0.013. The bulk 
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8S5 fit results are of good quality for both heating and 
cooling and yield a very small negative tte// that is not 
significantly different from fixing (Xeff to axY — —0.013 
as indicated by xt i2Dj- This indicates that the specific 
heat of the A^-Sm^ transition for the bulk LC behaves 
essentially as a clean XY transition consistent with the 
literature [23ll3ll |. However, it should be noted that the 
critical exponents for the parallel correlation length ) 
and for the smectic susceptibility (7) do not have XY 
values and that the A^-Sm^ transition for bulk 8S5 
exists in a complex crossover regime. 

As PS increases for the 8S5-|-sil samples, the fits for 
both heating and cooling are all of good quality even 
up to the highest density studied. In general, the ad- 
justable parameters have approximately a 10% uncer- 
tainty in magnitude for the bulk and the two lowest den- 
sity aerosil samples, except for Tc which has an uncer- 
tainty of approximately ±5 mK for all fits. This uncer- 
tainty grows for the higher density samples because an 
increasing temperature gap about the ACp peak of data 
is excluded from the fits due to rounding. These fits 
show that the critical character described by the expo- 
nent tte/ / of the A'^-Sm^ transition remains unchanged in 
the 8S5-|-sil samples for all ps, even above the ps = 0.1 
where previous studies of other LC-|-sil samples found no 
critical behavior 0, . It is apparent from the fits that 
the increasing temperature gap (itmin) is dominated by 
rounding on the low-temperature side of the SCp{NA) 
peak. This increasing temperature gap for the fits has 
the effect of causing the increase in the magnitudes of 
the correction-to-scaling coefficients although they 
remain well behaved. Fixing the ratio / D~ = 1 or 
setting = D~ — resulted in significantly poorer 
fits for all samples. It is also evident is that the critical 
amplitudes above and below the transition, A'^ system- 
atically decrease with increasing ps although the ratio 
A~^ /A~ remains unchanged. 

The scaling analysis for the A^-Sm^ transition in 
LC+aerosil systems follows that of reference |13| by us- 
ing the correlation length power-law to equate the cut-off 
length scale (maximum correlation length) S,m to a min- 
imum reduced temperature. This is closely related to a 
finite-size effect, however, there are no true pores and the 
LC filled voids are highly interconnected. The correlation 
length is given by 



(7) 



Since Eq. {7)) is defined only for T > T*, the finite- 
scale induced rounding of the transition is estimated in 
terms of the minimum reduced temperature on the high- 
temperature side of the transition t^, by the form 



6T*/T* « 24 



6w_ 



= 2 



2n 
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I/I'll 



(8) 

where the cut-off length scale is written in terms of n, 
the number of mean- void lengths Ig. The heat capacity 
maximum at the transition is given by substituting t+ 



into Eq. lO and is also explicitly defined for the T > T*. 
The explicit form for hm is 
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(9) 

Finally, the transition enthalpy SHna resulting from this 
scaling analysis is determined by using to truncate the 
integration of Eq. Q both above and below the transition 

Plotted on Fig.0]are the scaling trends for the A^-Smj4 
transition of 8S5-|-sil samples using the bulk A^-Sm^ 8S5 
critical parameters. Two choices for the cut-off length 
scale are shown. The first choice for the cutoff corre- 
lation length uses the mean distance between silica 
surfaces (mean void size) Iq = 2/aps, where a is the spe- 
cific surface area Q ■ The second choice allows S,m to vary 
as some multiple of lo, i.e. = nlo- The results for 
of the 8S5-|-sil samples are in very good agreement using 
= 3Zo, which more closely matches the measured sat- 
urated smectic correlation length However, ST*/T* 
appears sharper and SHna smaller for ps > 0.1 than 
predicted by this analysis. 



B. The Landau mean-field Smyl-SmC scaling 
analysis 

The excess specific heat due to the Smvl-SmC phase 
transition is shown in Fig. O for bulk 8S5 and 8S5-f sil 
samples over a temperature range down to Tac — 4.5 K. 
The bulk Smj4-SmC transition is well described by an 
extended Landau theory [2^ [s^ given by 

for T > Te 

(10) 

^ < 

where A is the dCp{AC) maximum at the transition Tc 
and r„i is the upper stability limit of the SmC phase. 
Results from fitting Eq. H10|) to the excess specific heat 
of the Sm^-SmC phase transition SCp for bulk 8S5 and 
8S5-t-sil samples over the temperatures T < Tac are tab- 
ulated in Table Hvl These results are the average of heat- 
ing and cooling scan fits. The stability limit of the SmC 
phase appears to shift upwards from the bulk value by 
approximately 1.6 K with increasing ps- Although in- 
triguing and indicating a stabilization of the SmC phase 
with increased ps, given the aging of this particular LC, 
it is not clear how robust is this result. The fit coeffi- 
cient A of Eq. pO|l appears to smoothly decrease with 
increasing ps- This is consistent with the behavior of the 
Sm^-SmC transition of 70.4 in hydrophobic aerosil [l3l| . 
The parameter A is also a measure of the tricritical na- 
ture of the Sm^-SmC transition. A log- log plot shown 
in Fig.|Slof A versus ps reveals a power-law scaling given 
by ^ « 
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TABLE II: Heating scan summary of the results of fitting Eq. Q to the excess specific heat peak ACp of the A''-Smyl phase 
transition on SSS+sil samples. The transition temperature Tc is given in kelvins, while the parameters Be and axe given 
in J K"^ g-'^. The parameters are dimensionless. All scans were fit from tmax = 10"^ to zttmin- AH parameters were free 
to vary in the fit except when the exponent was fixed to q?xy = —0.013 (denoted by the square brackets). 



PS 


Tc 




Be 


A+ 


A- 


D+ 


D- 


tmin X 10 


xi 





335.220 


-7x10-^^ 


1023.600 


-1023.760 


-1023.630 


0.0002 


0.0022 


+1.80/-12.3 


1.266 




335.222 


[-0.013] 
-1x10-* 


6.053 


-6.260 


-6.122 


0.0026 


0.3712 




1.319 


0.025 


334.712 


536.893 


-536.971 


-536.845 


0.0007 


0.0045 


-f3.50/-29.9 


1.073 




334.711 


[-0.013] 
-1x10"* 


4.757 


-4.880 


-4.741 


0.0450 


0.5234 




1.027 


0.041 


335.111 


383.615 


-383.643 


-383.536 


0.0016 


0.0060 


-I-2.56/-14.2 


1.287 




335.110 


[-0.013] 
-2x10"* 


4.124 


-4.194 


-4.076 


0.1117 


0.5715 




1.307 


0.051 


334.938 


154.089 


-154.083 


-154.025 


0.0045 


0.0086 


+ 11.4/-49.2 


1.134 




334.937 


[-0.013] 
-IxlQ-* 


3.228 


-3.249 


-3.186 


0.1954 


0.4153 




1.137 


0.078 


335.123 


376.628 


-376.688 


-376.614 


0.0010 


0.0036 


+3.08/-14.9 


1.163 




335.123 


[-0.013] 
-1x10"* 


4.049 


-4.146 


-4.064 


0.0628 


0.3348 




1.183 


0.105 


335.297 


259.368 


-259.267 


-259.171 


0.0049 


0.0095 


-K3.63/-19.3 


1.076 




335.297 


[-0.013] 
-7x10"* 


2.767 


-2.698 


-2.593 


0.4513 


0.9716 




1.084 


0.161 


335.162 


18.744 


-18.608 


-18.580 


0.0655 


0.0615 


+5.67/-42.3 


1.036 




335.161 


[-0.013] 
-4x10"* 


1.335 


-1.212 


-1.181 


1.0162 


0.9862 




1.037 


0.220 


335.498 


10.314 


-10.099 


-10.048 


0.1651 


0.1875 


+3.59/-57.5 


1.004 




335.496 


[-0.013] 
-6x10"* 


0.630 


-0.428 


-0.373 


3.9745 


5.2388 




1.004 


0.347 


335.145 


6.892 


-6.724 


-6.703 


0.1839 


0.1400 


+0.27/-85.1 


1.056 




335.144 


[-0.013] 


0.573 


-0.414 


-0.391 


3.0646 


2.4651 




1.056 



TABLE III: Cooling scan summary of the results of fitting Eq. @ to the excess specific heat peak ACp of the A'^-Sm^ phase 
transition on 8S5+sil samples. The labels are the same as those used in Table ITTl as well as the uncertainties of the adjustable 
fit parameters. 



ps 


Tc 




Be 


A+ 


A- 


D+ 


D- 


imin X 10 


xl 





335.198 


-7x10-^ 


1047.600 


-1047.800 


-1047.670 


-0.0001 


0.0019 


+3.42/-9.91 


1.463 




335.197 


[-0.013] 
-1x10-* 


6.565 


-6.826 


-6.688 


-0.0499 


0.2870 




1.532 


0.025 


334.632 


474.449 


-474.466 


-474.327 


0.0012 


0.0060 


-H7.86/-32.8 


1.159 




334.630 


[-0.013] 
-3x10"* 


4.422 


-4.484 


-4.334 


0.0956 


0.6721 




1.176 


0.041 


335.108 


141.506 


-141.494 


-141.395 


0.0056 


0.0167 


-h4.37/-23.5 


1.203 




335.108 


[-0.013] 
-1x10"* 


3.844 


-3.868 


-3.760 


0.1782 


0.6412 




1.207 


0.051 


334.845 


216.789 


-216.692 


-216.619 


0.0047 


0.0091 


-I-11.0/-36.4 


1.134 




334.843 


[-0.013] 
-2x10"* 


2.611 


-2.541 


-2.462 


0.3824 


0.8130 




1.138 


0.078 


335.057 


195.069 


-195.050 


-195.018 


0.0050 


0.0069 


+32.6/-32.7 


1.188 




335.055 


[-0.013] 


3.800 


-3.809 


-3.775 


0.2461 


0.3545 




1.192 


0.105 


335.179 


-0.0043 


11.848 


-11.888 


-11.823 


0.0425 


0.1152 


+37.0/-68.2 


1.145 




335.178 


[-0.013] 
-4x10"* 


4.239 


-4.301 


-4.232 


0.1061 


0.3208 




1.145 


0.161 


335.123 


48.578 


-48.455 


-48.427 


0.0192 


0.0211 


-I-16.4/-39.1 


1.151 




335.121 


[-0.013] 
-3x10"* 


1.647 


-1.538 


-1.508 


0.6025 


0.6853 




1.152 


0.220 


335.458 


33.646 


-33.468 


-33.415 


0.0363 


0.0481 


-H9.7/-99.6 


1.032 




335.455 


[-0.013] 
-3x10"* 


1.205 


-1.043 


-0.985 


1.1714 


1.6754 




1.032 


0.347 


335.14 


9.953 


-9.734 


-9.706 


0.1367 


0.1229 


-h5.65/-74.2 


1.129 




335.14 


[-0.013] 


0.332 


-0.119 


-0.088 


11.6406 


14.1668 




1.129 



IV. DISCUSSION 

Results have been presented from a series of high- 
resolution ac-calorimetric experiments on 8S5-|-sil disper- 
sions through the A^-Sm^ and Smj4-SmC phase transi- 
tions as a function of aerosil density. In bulk 8S5, the A^- 



SmA is a nearly clean XY-like transition with an a < 
while the Smj4-SmC transition is mean-field with a Lan- 
dau tricritical character. Our bulk measurements are 
fully consistent with these behaviors. The introduction 
of QRD in the 8S5-f sil system allows for the isolation of 
random-field, finite-size-like, and elastic strain effects at 
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TABLE IV: Summary of the results of fitting Eq. 1101 to the excess specific heat peak SCp of the SmA-SmC phase transition 
for 8S5+siI samples. The fit results have been averaged between heating and cooling scans and were done from ~ — 7 K below 
Tac to a point slightly below the rounded 5Cp peak. The transition temperature Tc and the SmC stability limit Tm are given 
in kelvins while the coefficient A is given in J g~^. No uncertainty is quoted for Tc as it was fixed to this final value for 
the last fit iteration. 



Ps 


Tc 


T 


A 


xl 





328.957 


329.26 ± 0.04 


0.108 ± 0.001 


1.187 


0.025 


328.706 


330.32 ± 0.38 


0.075 ± 0.002 


1.208 


0.041 


329.090 


330.61 ± 0.36 


0.056 ± 0.002 


1.328 


0.051 


329.010 


332.23 ± 0.31 


0.049 ± 0.003 


1.323 


0.078 


329.495 


330.49 ± 0.25 


0.051 ± 0.002 


1.258 


0.105 


329.250 


331.08 ± 0.40 


0.050 ± 0.010 


2.317 


0.161 


329.390 


329.97 ± 0.44 


0.025 ± 0.002 


1.160 


0.220 


329.860 


331.41 ± 0.19 


0.023 ± 0.003 


1.206 


0.347 


329.414 


330.37 ± 0.47 


0.016 ± 0.002 


1.211 



a pure XY and mean-field transition. 

The iV-Sm^ transition transition temperature does 
not exhibit a sensitivity to the aerosil (and in fact in- 
creases slightly) in stark contrast to that seen in nearly 
all other LC-|-sil studies. Although the ACp(iVA) de- 
creases uniformly in size with increasing p5, it remains 
sharp and is well characterized over the whole range of 
aerosil concentration. The power-law fits reveal that 
quasi-critical behavior is preserved and that the expo- 
nent tte// remains slightly negative and constant for all 
PS studied. See Fig. The systematic decrease of the 
coefficients but nearly constant ratio A'^ /A^ is a re- 
flection of the uniform decrease in the size of ACp {NA) 
(as well as 5Hma)- This analysis is consistent with the 
view that the aerosil gel effectively increases the nematic 
range. That is the nematic susceptibility decreases with 
increasing ps and this is reflected in a substantial change 
in the pseudo-critical properties for cyanobiphenyl-aerosil 
systems. For 8S5-f sil where the bulk material already has 
a long nematic range, little change is observed. 

Because of the good quality power-law fits available 
in this study, a detailed scaling analysis was performed 
and compared to the heat capacity maximum hm — 
AC™°'^{NA), transition enthalpy SH^a, and transition 
temperature rounding 6T* /T* of the A'^-Smj4 phase tran- 
sition. This analysis uses the bulk parameters and an ad- 
justable cutoff length scale S,m and found excellent agree- 
ment for h,n using ^A/(8S5-|-sil) = 3lo but this does not 
completely describe the behavior of SHna nor 6T* /T*, 
especially for ps > 0.1. Interestingly, when this scal- 
ing analysis is applied to the 8CB-|-sil and 80CB-|-sil 
systems, equally good modelling of hm is obtained for 
CM(8CB-|-sil) = lo and ^M(80CB-f-sil) = l.blo- It is sur- 
prising that such excellent agreement is obtained over the 
whole range of ps explored for these three LCs despite the 
apparent violation of the classic expectation of finite-size 
scaling (that is, the truncation of the bulk ACp behavior 
at the cut-off length-scale) . Also, the smectic correlation 
lengths measured for these systems are much larger than 
this cutoff length scale [s^ • This makes the geometric in- 
terpretation of the cut-off length puzzling but ultimately 



connected to the strength of the disorder. 

For the mean-field Sm^-SmC phase transition, the 
transition temperature also exhibits very weak sensitivity 
to the presence of aerosil (even increasing slightly rather 
than decreasing). This is also in stark contrast to that 
seen for the Sm^-SmC transition in 70.4-|-sil and 8S5 
in aerogel as well as the Sm^-SmC* transition of 
CES-fsil. The heat capacity associated with the Sm^- 
SmC transition in 8S5-|-sil presented in this work exhibit 
a systematic evolution from a Landau tricritical peak to a 
simple mean-field step with very little smearing observed 
for T > Tac, and then only at the highest ps- Good qual- 
ity fits were made using an extended Landau form that 
found the upper stability temperature increasing with 
Ps consistent with the shift of the observed transition. 
In addition, the coefficient A, representing the heat ca- 
pacity maximum at the Sm^-SmC transition, exhibits a 
scaling with the aerosil conjugate density as A oc Ps'^'^, 
see Fig. |S1 It is clear that this transition, while evolv- 
ing in a systematic fashion, remains mean-field over the 
entire range of ps studied. This may be a reflection of 
the Sm^-SmC transition behaving as if it were at its up- 
per critical dimension moderating the effect of the QRD 
produced by the aerosil gel. 

In this study on 8S5-|-sil, it appears that the random- 
field QRD of the aerosil gel plays a relatively weak role 
at the A^-Sm^ transition (due to the non-polar molecules 
and the long nematic range) nor at the Sm^-SmC tran- 
sition (due to either its proximity to its upper critical 
dimension or that the QRD does not couple linearly to 
the SmC order parameter and so, does not provide a 
random-field interaction) . It is clear that a special length 
scale is present with > lo but it does not appear to 
play a leading role for the Smvl-SmC transition. How, 
then, does one understand the evolution of this transition 
with increasing QRD of aerosil and in light of the other 
LC+sil and LC-|-aerogel results? One possibility is that 
the aerosil gel, due to its fiexibility, is closer to thermo- 
dynamic equilibrium with the host LC and so causes the 
LC-f sil to behave as a stiffer LC as well as provide QRD. 
The aerosil gel has been shown to exhibit dynamics cou- 
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FIG. 4: Scaling analysis of the 7V-Smj4 phase transition in 
8S5+sil (solid circles) on log-log scales. Top Panel: Tran- 
sition enthalpy SHna- Middle Panel: Excess specific heat 
maximum hm = AC]^^^. Bottom Panel: Transition round- 
ing in reduced temperature Stmin = \ + The bulk 
scaling predictions are given by Eqs. Q and JHJ where the 
mean-void length lo is used as the cut-off length scale (solid 
lines) or 3Zo (dashed lines). 



pled to the host liquid crystal |3J| and recent work has 
followed its quenching as aerosil density increases [35, 36]. 
This increase in the effective microscopic elastic stiffness 
of the LC would be similar to the engineered stiffening of 
polymer composite materials. A consequence of the stiff- 
ening of the nematic phase in LC-|-sil systems has already 
been discussed previously P, Q, 0, Is^] in terms of the 
decrease in the nematic susceptibility with increasing ps 
to explain the crossover from tricritical to XY behavior 
for the A'^-Snij4 transition. The present view is to extend 
this concept to the general stiffening of the LC within 
the LC+sil. 

In the case of the Sm^-SmC transition the quenched 
disorder is also changing the character of the transition. 
With increasing aerosil density, the peak at the transi- 
tion is being suppressed in favor of a step-like behavior. 
This indicates that the transition is being driven away 
from tricriticality toward mean-field character alone. The 
proximity to tricriticality is controlled by the layer com- 
pression elastic constant, B, in the modified fourth-order 



P, (gem" ) 



FIG. 5: Scaling plot of A (SCp at T = Tac), obtain from 
fitting the excess specific heat of the Smj4-SmC phase tran- 
sition to Eq. 1101 . against ps- The solid line has a slope of 
—0.5 while the dashed line a slope of —0.56. 



free energy coefficient b ^ bo — 2X^/B, where bo is the 
original value of this coefficient and A is the strength of 
the coupling between the tilt angle and the layer com- 
pression 0. The mean- field phase transition behavior 
and the decreasing size of the step in the heat capacity is 
indicative of a hardening of the layer compression elastic 
constant B. 

This effect would serve to explain the evolution of the 
Smj4-SmC transition from Landau tricritical to step in 
terms of a stiffening of the smectic layer compression 
modulus B and account for the general stability of the 
transition temperatures. Note that for 8S5, the nematic 
and smectic phases are almost fully decoupled in the 
bulk while there remains significant coupling in the 8CB 
and 80CB LCs. It is the latter two LCs that exhibit 
the strong non-monotonic transition temperature shifts 
downward with ps and so may be reflecting QRD effects 
at the I-N transition js^l- It would be expected that 
twisted (or chiral) phases would be more disordered and 
smeared due to this general LC stiffening as seen in a re- 
cent Smj4-SmC* in aerosil study 0. Also, if the silica gel 
is too rigid, then this LC stiffening effect is supplanted 
by the strong quenched disorder as seen in LC in aerogel 
and in high density aerosil gel. 

Finally, there seems to be an apparent inconsiste ncy 
between the experimental results of the x-ray study pUJ l 
and this calorimetry study on 8S5-f sil. In the x-ray work, 
it is clear that the temperature dependence of the tilt an- 
gle of the SmC phase is insensitive to ps, strongly sug- 
gesting that the order-parameter tp for the SmC phase 
remains unchanged from bulk for all 8S5-|-sil samples. 
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FIG. 6: Comparison of effective specific heat critical expo- 
nent Qe// obtained from fitting Eq. Q to the A/'-Smj4 phase 
transition data for 8S5+sil (sohd circles) and SCB+aerosil 
(open circles) samples. 



However, if the system remains mean-field, then the re- 
lationship Cp oc di/j^/dT should hold. In this case, the 
dramatic variation of SCp{AC) with ps suggests that the 
SniC order-parameter should be strongly ps dependent. 
Comparison of the 5Cp data in Fig. |21 (one data point 



every 0.1 K) and the tilt angle data in Fig. 6 in ref. 
(one data point every w 1 K) suggest that the situation 
is more ambiguous. A suppression of the tricritical char- 
acter may be consistent with the more sparse x-ray data 
(see note [21] in ref. [lol|). 

Clearly, there is an important need for theoretical ef- 
forts at understanding quenched random disorder on the 
Sm^-SmC transition. The beginnings of such a the- 
ory, although for the anisotropic aerogel disordered Sm^- 
SmC transition 'ssj, are emerging but much is left to be 
done. A coherent framework must account for the appar- 
ent connection with elasticity of both the host material 
and the disorder inducing gel. It must also provide for 
an interpretation of the cut-off length scale that seems to 
play an important role for the A^-Sm^ as well as the ob- 
served scaling of the ACp (AC) peak with ps ■ Additional 
x-ray and calorimetric experimental work with anisotrop- 
ically aligned Sm^-SmC as well as Sm^-SmC* in aerosil 
gel systems are also needed and would provide important 
data as to the role the polar nature of the LC and the 
local aerosil-LC interaction. 
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